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Abstract. The double and single differential cross-sections with respect to positron 
and electron energies as well as the total cross-section of triplet production in the 
laboratory frame are calculated numerically in order to develop a Monte Carlo code 
for modelling electron-photon cascades in a soft photon field. To avoid numerical 
integration irregularities of the integrands, which are inherent to problems of this 
type, we have used suitable substitutions in combination with a modern powerful 
program code (Mathematica) allowing one to achieve reliable higher-precission results. 
The results obtained for the total cross-section closely agree with others estimated 
analytically or by a different numerical approach. The results for the double and 
single differential cross-sections turn out to be somewhat different from some reported 
recently. The mean energy of the produced particles, as a function of the characteristic 
collisional parameter (the electron rest frame photon energy), is calculated and 
approximated by an analytical expression that revises other known approximations 
over a wide range of values of the argument. The primary-electron energy loss rate due 
to triplet pair production is shown to prevail over the inverse Compton scattering loss 
rate at several ('~-"2) orders of magnitude higher interaction energy than that predicted 
formerly. 



1. Introduction 

There are two main reasons wliy triplet pair production (TPP) has been commonly 
ignored in astrophysical apphcations. The first reason is that TPP is a third-order 
QED process. The second reason is the extremely complicated and long expression for 
the total differential cross-section [|1|. This considerably complicates the modelling of 
the energies and momenta of the final three particles in comparison with the case of 
considering only the two major processes for electron-photon cascade in a radiation field: 
pair production and inverse Compton scattering. Apart from the formal complications 
there are serious mathematical problems to be overcome connected with numerical 
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calculation of the double differential cross-section (DDCS) and single differential cross- 
section (SDCS) in the laboratory frame. A typical problem is the integration over the 
cosine (cos^_) of the polar angle ^_ of the produced electron momentum p_, where both 
the integrand irregularities coincide with the integration limits whose semi-vicinities 
provide the major contribution to the integral. 

Despite the above-mentioned arguments, at ultrarelativistic electron energies TPP 
becomes a prevailing process playing an important part in the electron-photon cascades 
in a soft background photon field that form the energy spectrum from a variety 
of astrophysical sources. This fact has been recently emphasized and confirmed by 
incorporating TPP into full cascade calculations HJ-p. 

Based on the recently revived interest in a more precise simulation of electron- 
photon cascades in a photon field, we have started to develop a Monte Carlo code 
for modelling TPP. Our intention is to use this code for a more detailed study of the 
development of electromagnetic cascades in thermal fields. To realize this intention we 
have decided to follow a method like that suggested in paper 0]. Thus, as an initial 
step we have precisely recalculated the DDCS and SDCS with respect to electron and 
positron energies as well as the total cross-section of TPP in the laboratory frame. An 
obstacle to overcome here is the presence of the above-mentioned integrand irregularities 
in combination with extremely short integration intervals. So, one purpose of the present 
study is to search for ways to avoid these intrinsic difficulties and to achieve more precise 
results. Another purpose of the study is to calculate and analytically approximate 
the mean energy of particles produced as a function of the characteristic collisional 
parameter. Investigating the primary-electron TPP energy loss rate is also an important 
aim of the work. 



2. Method and results 



2.1. Calculation approach 

Let us first consider the basic expressions of interest corrected for typographical 
errors that have appeared in many papers. The DDCS with respect to the positron and 
electron energies is given by 



dV 

-{Eq, eq, s, E+, EJ) 



-p+p- dx dy[ai{yn,e.x - y){y -Vmin)] ' \ Xd0+, (1) 

^min ?/min 

(x) Jo 



where p+, Qj^ and p_, are momenta, energies and polar angles of the produced 

positron and electron, respectively, x = cos 9^ and y = cos9_, 0+ is the azimuthal 
angle of the positron, Eq and Eo are the energies of the incoming electron and photon 
respectively, s = EqEo{1 — Pcos9) is the characteristic collisional parameter {9 is the 
collision angle) representing the photon energy in the electron rest frame (ERF), and X 
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Figure 1. Integrand function ^(iJo, eo, s, S+, x, ?/, 0+)[ai(2/niax-2/)(2/-ymin)] ^'^ 
in arbitrary units (before change of variables) versus the polar-angle cosine y = cos 6'_ 
at Eo = 10^ £0 = 10-^ s = 10^ E+ = 4.999925000425 x 10^, E- = 2.500037499837 x 
10^ X = cos 9+ = 0.999999999995, 0+ = n. 

is a cumbersome expression that is given in the appendix; the quantities «/ = e^/{hc) 
and To = e^/inieC^) are the fine structure constant and the classical electron radius, 
respectively, and irie is the electron rest mass. The limits of integration are: 

1/max = [h + {hi - aiCi)^/^]/ai, 

?/min = [h - {h\ - aiCi)^/^]/ai, 

Xrnln = -{F^ + P-FI'^) / {p+Ptot) . 

where 

oi = pt {Ptot + P+ - 2PtotP+x) , 6i = Ap_ (p+x 

a2 2 2 • 2 /I 

Ci = ^4 — sm 6'+, 

A = l + s + p+PtotX - Etot{E+ + + E+E_, 

Fi = E'i — {Etot — E+)E^ + s — EtotE+, 

F2 = {E+ + E_ — Etotf' — 1, Etot = Eq + Sq, 

"jf and k are the momenta of the incoming electron and photon, respectively. Let us 
note that throughout this paper the energy quantities are in units mgC^. Note that 
the integrand in equation (|I]) has two irregular points with respect to the variable y. 
These points, y = ?/min and y = ^max, coincide with the integration limits. Also, the 
integration intervals over x and y, especially at high energies, are extremely short. In 
addition, even within such narrow integration intervals the integrand changes sharply 
with changing of y (see figure 1). Because of these peculiarities of the integrand, an 



Ptot), 



and 



'^tot — 1^0 + k 
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Figure 2. Integrand function 2X{Eo,eo, s,E+, E-,x, rj, (f>+)/ {l + r/'^) in arbitrary units 
(after change of variable y) versus the new variable rj at the same fixed values of the 
remaining variables as in figure 1. 



accurate calculation (numerical integration) can be successful only when a sufficiently 
high-precision number is used. This was first been pointed out by Mastichiadis 
who underlined the necessity of quadratic precision to calculate DDCS and SDCS, 
and reconsidered his own formerly obtained results [0]. Since our purpose here is to 
perform similar calculations with a higher precision (e.g. up to 80 significant decimal 
digits) we have used the program code Mathematica [|^ , which allows one to work with 
arbitrary high-precision numbers. Thus, one can both eliminate the precision number 
conditioning and ensure a reliable precision of the final results. In addition, this code has 
an adaptive program for quadrature of multiple integrals that precisely approximates the 
fast-changing integrand and permits one to obtain results with a prescribed precision. 
Nevertheless, the extraordinary character of the integrand in equation (^ leads to a 
fast growth of the CPU time with the increase of the interaction energies. Besides, 
it is not unknown for the program to fail in some cases. As a result of searching for 
ways to eliminate the above-mentioned problems, we came to the following change of 
variables that led to acceptable integration intervals and acceptable smooth behaviour 
of the integrand: 



Equation (|25|) is, in fact, one of the Euler substitutions that is appropriate for this 
case and leads to the integral: 






{2b) 
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Figure 3. DDCS (a/r§) M^CT/(d£'+di?_) as a function of and £_ at various 
values of s; Sq = 10~^, Eq = s/eq (glancing collision). 




E+ E+ E+ 



Figure 4. Three-dimensional contour plot of {afrQ)~^(Pa/{dE+dE-) as a function 

of i?+ and E-, at values of s as in figure 3, respectively, around the first peak of figure 
3 disposed near i?+,inm and E'-^mm; eo = 10~^, Eq = s/eq (glancing collision). 

/■27r 

X / d(/.+X(£;o,£o,s,£^+,£^-,e,?7,0+)/(l + (3) 

dV/(d£;+d£;_) = (a/ ro7(47r2a;/'s)Kp_7(£;o,£o, 

It can be seen that the new variables lead to an enlargement of the integration scale 
and removal of the integrand irregularities. As one may expect, the integrand becomes 
a smooth function of 77 (figure 2), which leads to an acceleration of the calculation 
procedure, increase of the precision, and elimination of any computational failures when 
Mathematica is used. 
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Figure 5. Dependence of the quantity (afr^) ^D{s,z) = (£'_ — E^,niin){E+ — 
^^+,min)[dV/(d£'+d£'_)] on the parameter z ^ {E_ - E_^„,in)/{E_^„ied - E_^,^in) 
[E-^rned = 0.5(£'-,min + -E-anax)] for i?-,min < E- < E-^rned and various Outgoing 
positron energies E+ = W{a),lQ^{h),lQ^{c),l(f{d) and lO^(e); E^^ = 10^ and Eq = 
10-3 (s ^ 105). 

2.2. Double differential cross-section 

We have obtained precise results for DDCS as a function of £"0, ^o, s, i?+ and E^. Some 
of these results are shown in figures 3 and 4 in a three-dimensional surface form and a 
three-dimensional contour form, respectively. It is seen (figure 3) that the dependence 
of d^(j/(d-E'+di?_) on and considered over the whole range of the arguments 
and at various fixed values of the parameter s = EoEq ( in the case of a glancing 
collision when 6 = it /2 ), is represented by a double-peak surface whose secants with the 
planes E+ = const are symmetric curves with respect to the point {E+, E^^rned), where 
E-^med = 0.5(£'_^min + -E'-,max)- With the iucrcase of s, the surface peak heights increase, 
and the surface itself (as well as the corresponding dependence of DDCS on E^ and 
becomes sharper. Also, up to values of s = 10^, the peaks change their positions over the 
plane {E^,E_}. So, the common coordinate of both peaks along E^ axis is shifted to 
lower values of E+. In addition, the first peak (disposed below the point {E+, E^^med)) 
is shifted to lower values of E_ (see figure 4), and the second one, to higher values of 
E_. At values of s increasing above s = 10^ the peaks do not change their positions; the 
surface is as if consisting of two spikes (with a common E^ coordinate) whose positions 
are close to the points (ii^+^min, i?_^min) and (i^+^min, -E-,max), respectively. The same 
results as above are represented in figure 5, but in a parametrized form proposed in @], 
where the quantity D{s, z) = {E_ — E_^-^i^){E^ — i?4._mm) [d^cr/(d-E+d_E_)] is considered 
as a function of the parameter z = (-£/— — £'_^inin)/(ii'_^me(i — ,min 

) at fixed values of 

E^ and s. It is shown in that D can be considered as dependent only on s and z 
if £^0 ^ 1 ^ ^0- Such a parametrization is useful for application to a Monte Carlo 
code for modelling electron-photon cascades in a soft photon field taking into account 
the contribution of the TPP process. Then, on the basis of tabulated data D{s,z) 
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Figure 6. Dependence of the quantity (oifrl) ^C{s,y) = {E+ — i?+^niin)[dcr/di?+] 
on the parameter y = (Ej^. — i?+,min)/(-E'+,max ~ -E'+.min) for various ERF energies of 
collision s = I0{a),10^{b),l0^{c),10'^{d),l0^{e),10^if),10'^{g) and I0^{h). 



Table 1. TPP total 


cross-sections (in units af 
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s,s± 


^tot^Haug 


(^tot,Mast 


<^tot,Our 


^toti,Our 


4.01 


5.80 X lO-'^ 


5.4 X lO-'^ 


5.8 X lO-'^ 


0.092 


4.1 


7.53 X 10-5 


7.6 X 10-5 


7.5 X 10-5 


0.102 


5 


0.0170 


0.019 


0.0170 


0.221 


7 


0.179 


0.19 


0.179 


0.56 


10 


0.594 


0.59 


0.59 


1.12 


102 


7.21 


7.3 


7.3 


7.9 


10^ 


15.1 


15.3 


15.1 


15.8 


lO" 


22.6 


22.4 


22.7 


23.2 


10^ 


29.9 


29.7 


29.9 


30.5 


lO*' 


37.1 


37.0 


37. 


38. 


10^ 


44.2 




44. 


45. 


108 


51.4 




52. 


52. 



one can determine DDCS at any combination of Eq, Sq, s and positron energy Ej^ @]. 
Comparison between our results and those obtained in 0] shows that our curves pass 
through maximum and that at lower electron energies tending to i^.^min, where the 
calculation is sensitive to loss of precision, they essentially fall below the corresponding 
curves obtained in 0. 

2.3. Single differential cross-section 

The SDCS do"/di?+ has been calculated by integrating d^cr/(di?+d-E'_) over E^ with 
integration limits E_\^^ = Q.b[Et,t- E+±{Ptot-p+){l-2/ Bf'% {B = l + s-EtotE+ + 
PtotP+)- In the same way, dcr/dE_ has been calculated by integrating cPa/ (dE^dE^) 
over The use of an optimum-power spline technique allowed us to obtain precise 
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results for sufficiently high values of s > 10*^. (We consider as an optimum power of the 
spline that one, above which the results from the integration remain stable.) The same 
results have been obtained by direct integration (without spline interpolation) but using 
considerably longer CPU time. Some of the results obtained are shown in figure 6 in a 
parametrized form where the quantity C{s,y) = — _E+ min)[d(T/dii^+] is considered 
as a function of the parameter y = — -E+^min) / (-E+,max — -E+,min) at fixed values of 
s. For Eq ^ 1 ^ eq, C{s,y) depends only on s and y, and (when tabulated) allows 
one to determine the SDCS for any combination of Eq, Eq and s [Q. The comparison 
between our results and those obtained in shows that the corresponding curves have 
the same behaviour with a characteristic maximum. The difference is that the left- 
hand part of each curve of ours (including the maximum), where the calculation is more 
sensitive to loss of precision, is as if shifted right with respect to the analogous part of 
the corresponding curve obtained in [^. 

2.4- Total cross-section 

We have also calculated the total cross-section atot by integrating da/dE^ over E+. The 
integration limits are: E+\^f^ = {Etot{s -I) ± Ptot[s{s -4:)Y^^}/{l + 2s) §. The results 
obtained are compared, in table |l|, with the results obtained by other authors @], 
The agreement is excellent and may be considered as an indirect confirmation of the 
precise character of our calculations. 

2. 5. Mean energy of produced particles 

Finally, we have calculated the mean energy E^ = -E'+,m = of the produced 

particles on the basis of the relations: 



The results obtained for -E+^m and E^^m are practically coincident, which may be 
considered as another conffimation of the precise character of the calculations performed. 
On the basis of the parametrization approach developed in 0] we can show that 
the ratio E^/ Eq is a function only of s when Eq ^ 1 ^ eq. In this case, the 
integration limits -E+.min and -E+,max (see above) are expressible as -E+,min = Eofmin^s) 
and ii^+,max = EQfma.x{s), where the functions /min and /max depend only on s. After the 
change of variable E^ = xEq, instead of equation (H) we can write: 




(4) 



E, 



'm 



dx.x.[(d(j/dx)(i?o, ^0, s, xEq)]. 



(5) 



/min{s) 



In the same way, taking into account that Otot depends only on s, we obtain 
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Table 2. Normalized mean energies ol produced particles. 







EmiS±/Eo = Ejni£0 


0.195 ln2(2si) 
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U.oO 
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1.44 


1 Q A 

1.84 


1 OQ 

i.Uo 


7 


1.68 


2.10 


1.36 


10 


1.97 


2.39 


1.75 


102 


5.00 


5.50 


5.47 


10^ 


10.1 


10.7 


11.3 


10"* 


17.7 


18.5 


19.1 


10^ 


27.1 


28.2 


29.0 


lO*' 


39.4 


40.8 


41.0 


10^ 


54.2 


55.8 


55.1 


108 


70.0 


71.5 


71.2 



^/max(«) ^/max(s) 

CTtotis) = dx[{da/dx){Eo,eo,s,xEo)] = dx[(da-/dx)(Eo, 4, s, xEq)] (6) 

where i?Q and are some other values of the incoming electron and photon energies, 
respectively, but such that the value of s is retained. Based on equation (||) we may 
conclude that {da/dx){EQ, e'q, s, xE'q) = {da/dx){Eo, Eq, s, xEq), and consequently (see 
equation d)) E^/Eq = E'^/ E'q = f{s), where 

f{s) = arol{s)N{s), (7) 

and 

/•/max(s) 

N{s) = / dx.x.[{da/dx){Eo,eo,s,xEo)]. (8) 

-^/min(s) 

Thus, the knowledge of /(s) allows one to determine E^ for any pair of values of 
Eq and s. At fixed values of Eq and eo, /(s) describes, in practice the dependence of 
Em on the angle of collision 9. 

The results calculated for E^s/Eq versus s are represented in figure 7(a) by black 
squares (see also table H). At s > 10^ they are well described by the dependence (curve 



h{s) = Ern{Eo,s)s/Eo = 0.195 ln2(2s) 



(9) 
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Figure 7. (a) Plot of the quantity /i(s) = Ems/Eo (tlie mean energy Em of produced 
particles normalized to Eq/s) versus ERF collision energy s. The black squares 
represent the results from our calculations fitted by a cubic spline (curve(6o)); at s > 
10^ curve (bo) is well described by the function /i(s) — 0.195 ln^(2s). Curves (eo) and 
(flo) correspond to the approximation /i(s) = 0.71s^/^ and /i(s) — \'a'{s/A)/atotBH 
based on the results of Dermer and Schlickeiser [9], and Feenberg and Primakoff 
[10], respectively, (b) Plot of the quantity Lpi{si_) — EmiSA^/EQ = Emi^o (the 
mean energy of produced particles in an isotropic and monochromatic soft photon 
field, -Emi, normalized to Eo/sj_) versus s± = SqEo- The black circles represent the 
results from our calculations. At s± > 10^ they are fitted by the function (pi{s±) = 
0.195 ln^(2sj_) represented by curve (&i). Curves (ei), (di), (ci) and (oi) correspond 
to the approximations ipi{s±) ~ f^l^^j V'llsj-) — 0.57s'];^'^, (/?i(sj^) — 2.5s*]^^^, and 
= 'PFp{s±) based on the results of Dermer and Schlickeiser [9], Mastichiadis 
et al [2], Mastichiadis et al [4] and Feenberg and Primakoff [10], respectively. 



The concrete calculations are performed at = 10^^, 6 = 7r/2, and various values 
of Eq leading to various values of s = EqEq. Nevertheless, the results obtained for 
Em/Eo should, as a whole, depend only on s independently of the concrete values of 
£■0, Eo and 6. Thus, on the basis of a special case we obtain the dependence Em{EQ, s) 
having a more general validity. In figure 7(a) we have also graphically represented two 
other estimates of the function /i(s) = Em^Eo, s)s/Eq obtained by other authors. The 
line (eo) corresponds to the estimate Em{EQ,s) = CTl-Eo-^'^'^obtained analytically by 
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Dermer and Schlickeiser |p. At values of s < 10^ this line passes through our points 
(squares). At values of s > 10^ the line goes far above our points, thus predicting several 
orders of magnitude higher results for E^. The curve (ao) corresponds to the estimate 
Em{Eo, s) = (Eo/s) It? {s / A) / atotBH{s) obtained on the basis of the theoretical approach 
developed by Feenberg and Primakoff |T^ (see equations (19), (21), (22) and (31) in 



(^totBH{s) is an analytical approximation to the Bethe-Heitler formula for atot{s) 
normalized to the quantity afr^. It is seen that curve (ao) lies essentially below our 
points and predicts several times lower results for Em{Eo, s). A reason for this is that 
curve (ao) describes an approximation obtained analytically as a lower limit of the true 
dependence fi{s). Another reason, that was pointed out by Mastichiadis et al ^ is the 
neglect (in |jlO|) of the recoil of the primary electron in the electron rest frame. 

In order to determine the mean energy Emi of a particle produced by relativistic 
electron-photon collision in an isotropic and monochromatic soft photon field one should 
additionally average E^ over the angle of collision 9. Then the expression of -£/rnj is 
obtained in the form: 

^™(^o,s±) =^o<^(s±), (10) 
where (p{s±) = Ni{s±)/atoti{s±) is a function of s± = EoEq; 

psoEoil+IB) 

N,{s^) = {2(3elE',y' / sN{s)ds, (11) 



4 



reoEoil+fS) 

atoti{sA_) = i2l3elEl)-^ j satot{s)ds. (12) 

The results calculated for ipi{s_[_) = Emi{EQ, s±)s±/ Eq = Ejni{Eo, s±)eo versus s± 
are represented in figure 7(b) by black circles (see also table At > 10^ they are 
fitted by the dependence (curve (6i)) 



<^i(s±) = Emi{Eo,s^)eo = 0.195 W(2sx) (13) 



(i.e. Emi{EQ, s±) Erai{£o,s±) = (0.195/£:o) ln^(2s^)) that has the same form as 
the dependence /i(s) (equation (||)). There are four more curves represented in 
figure 7(b), which describe some approximations of the function fi{s±) obtained by 
different authors. The line (ei) corresponds to the estimate ^i{s±) = fs^^ 0. 
Certainly, it almost coincides with the line (cq) in figure 7(a), and at s < 10 passes 



Numerical analysis of triplet pair production cross-sections 



12 



in immediate proximity to our points (circles). With the increase of s± above 10^, 
the discrepancy with our results also increases, achieving two orders of magnitude at 
s = 10*. The lines (di) and (ci) correspond to the approximations (pi{s±) = O.STs'}'^^ and 
^iis±) = 2.5s^^^ obtained by Mastichiadis et al in 1986 and 1991 0, respectively. 
The first approximation (line (di)) has been obtained by numerical calculations. It 
is near that obtained by Dermer and Schlickeiser, and has a similar behaviour with 
respect to our results. The latter approximation is obtained after reconsidering the 
first one and performing improved calculations and computer simulations. In the 
interval from s = 10^ to s = 10* the line (ci) lies just above our results. For 
completeness, we shall also briefly consider the results for (pi{sj_) = (Pfp{sj_) = 
(l/s_L){s_L[(lns_L - ln2 - 1)^ + 1] - [(ln2 + 1)^ + l]}/crtotiFP obtained on the basis 
of the approach of Feenberg and Primakoff (curve (cti)) by using the above-mentioned 
equations (19), (21), (22) and (31) in |10|; crtoUFP is derived from equations (19) and 
(21) in |jTO|, and is normalized to (divided by) a/r-Q. The corresponding curve ((ai) in 
figure 7(b)) almost coincides with the curve (oq) in figure 7(a), thus showing the same 
behaviour with respect to our results. The reasons for such a behaviour are pointed 
out above. The resemblance between our results and those of Feenberg and Primakoff 
is that V5i(s_l) (/i(s)) is obtained to be proportional to In^ s± and In s± (In^ s and In s) 
respectively, and not to a power function of s± (s) as in the other known approximations. 
Let us finally note that the lowest estimate of Em ^ 2/6o was mentioned by Blumenthal 



11 1 . Thus, it appears that there is a natural tendency to improve the calculation 



accuracy, leading to the results obtained here. 



2.6. Primary- electron energy losses 



The primary electron energy loss rate L^pp (the energy lost per unit time) due to TPP 
in an isotropic and monochromatic soft photon field is given by the expression 

Ltpp = 2cnEoNi{s±) = 2cnEmiiEo, s±)atoti{s±), (14) 



where c is the speed of light and n is the number of photons per unit volume. According 
to the results given in table |I]the values of crtoti{s±) at s± > 10^ are described correctly 
by the Haug formula 

28 218 

(rtoti{s±) = afrl[—\n{2s±) - —]. (15) 



Also, as shown above, at sj_ > 10^ the values of Emi{EQ, s±) or Emi{£o, s±) are described 
correctly by equation (|T3|) . Therefore, based on equations (|13D-(|1^) we can write the 



following analytical expression of Ltpp normalized to the quantity x = cnnrQ/sQ: 
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Figure 8. Plot of the quantity qtpp = Ltpp/x (TPP energy loss rate Ltpp 
normalized to x — cnirr'^) versus s±_ = EqEq, compared with the analogous quantity 
Qics = Ljcsix concerning ICS energy loss rate Lies (curve (/)). Curve (&) is 
obtained on the basis of precise calculations performed in this work. Curves (c), (d), 
(e) and (a) represent some approximations based on the results of Mastichiadis et al 
[4], Mastichiadis et al [2], Dermer and Schlickeiser [9] and Feenberg and Primakoff [10], 
respectively. 



qTPp{s^) = Ltpp/x = 0.386a/ ln2(2sx)[ln(2sx) - — ]• (16) 

o4 

The primary-electron energy loss rate Ljcs due to inverse Compton scattering (ICS) 
in an isotropic and monochromatic soft photon field (normalized again to x) is given by 



g/C5(s±) = L/c5/x = ln(4s_L) - — . (17) 

6 

The normalized losses qrpp and qics versus sx_ are compared graphically in figure 8. 
There, curve (/) represents the dependence qics{.s_i) given by equation (^). Curve (6) 
is obtained on the basis of precise numerical calculations performed in this work. It can 
be seen that, according to our results, TPP losses become prevalent and increase above 
ICS losses at values of s_l exceeding a threshold sth ~ 10^. Certainly, this threshold 
is considerably higher (five orders of magnitude) than another threshold s* = 250 at 
which the interaction lengths of TPP and ICS become equal Q]. An estimate of Sth 
derived from the equality qTPp{s±) = qics{s±) (by using equations (|T6|) and(p^) is 
Sth ^ 1-6 X 10^ 

The estimate of qTPp{s±) = (5/7r)Q;/sy^[^ ln(2s_L) — ^] obtained by Mastichiadis 
et al in 1991, is shown by curve (c). The threshold predicted in this case is Sth — 
10^. Curve (e) represents the estimate qTPp{s±) = ^a/s^^obtained by Dermer and 



Numerical analysis of triplet pair production cross-sections 



14 



Schlickeiser 0. It predicts a threshold sth — '2 x 10^. Two more estimates of 
the dependence qTPp{s±) obtained from the resuhs of [jTO| and are illustrated by 
curves (a) and (d), and give unrealistically high and low thresholds, respectively. So, 
the consideration performed here, of the results obtained by different authors for the 
primary-electron energy loss rate due to TPP, confirms the existence of a tendency to a 
permanent improvement of the calculation accuracy. Because of the efforts to overcome 
the calculation difficulties and additionally increase the calculation precision, one can 
accept the results obtained here as reliable and accurate. They show that the electron 
energy losses due to TPP, e.g. in cascading processes occurring in pulsars (|^, |l3l) or 
in photon background field (0,0], IHI), are lower than those formerly predicted. The 
differences between the results for DDCS, SDCS, Em{EQ, s),a.nd Emi{EQ, s±) obtained 
here and those obtained in other works might lead to differences in the results from 
modelling electron-photon cascading in a soft photon field. Certainly, the detailed 
simulations now in progress will reveal the influence on the final results of the differences 
and factors discussed in this paper. 



3. Conclusion 



In order to develop a Monte Carlo code for three-dimensional modelling of TPP, we have 
undertaken a series of systematic precise calculations of DDCS and SDCS with respect 
to the produced electron and positron energies, as well as of the total cross-sections of 
TPP in the laboratory frame. The behaviour of the mean produced-particle energies has 
also been investigated in detail. To avoid crucial irregularities and sharp variations of 
the integrand, and extremely short integration intervals in the expressions of DDCS, 
SDCS, total cross-section, and mean produced-particle energy, we have used some 
appropriate mathematical approaches such as suitable changes of variables, optimum- 
power spline technique etc. These approaches lead to simpler and regular expressions of 
the cross-sections as well as to stable, accurate and accelerated calculation procedures. 
In addition, the use of the modern powerful program code Mathematica, working with 
arbitrary precision numbers, allowed us to obtain reliable high-precision results. Thus, 
the DDCS, SDCS and total cross-section have been computed for a variety of initial 
and final parameters characterizing TPP. The results for the total cross-section are in 
excellent agreement with ones obtained by other authors. However, there are some 
discrepancies in the results for DDCS and SDCS that might lead to differences in the 
results from modelling. 

The mean produced-particle energy E^ is analytically confirmed (in a general 
form) to be proportional to the incoming electron energy Eq and to a function f{s) 
of the coUisional parameter s only, i.e. E^ = Em{Eo,s) = Eof{s). It is also 
confirmed analytically that in an isotropic and monochromatic soft photon field the mean 
produced-particle energy Emi (averaged over the angle of collision 9) is proportional to 
Eq and a function (p{s±) of the product s± = EqEq, i.e. E^i = Emi{EQ, s±) = Eq{p{s±). 
Such a general behaviour established of Em and Emi is in agreement with the results 
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of other authors (0, 0], 0, [[T^])- However, there are some essential differences that 
would also lead to different results from modelling. So, the mean produced-particle 
energy Em or Emi obtained here is proportional (apart from Eq/s or Eo/s±) to ln^(2s) 
or ln^(2s^), respectively (see equations (^ and ([T3|)). At the same time, some earlier 
investigations of this question have led to a proportionality to a power function such as 
s^^^ 10, s^jf'^ or s*};^^ [Q. The indicated differences in the determination of E^i lead 
to differences in the determination of that threshold level of s± = Sth above which the 
primary-electron energy losses due to TPP become prevalent over the ICS energy losses. 
It is shown here that the value of Sth ~ 10^, that differs from the values of Sth ~ 10^ 
0] or 2 X 10^ obtained formerly. The last-mentioned result means that there has 
been some overestimation of the role of TPP energy losses in some astrophysical studies 

(i-i, i, n, ii). 



Appendix. Expression of the integrand function X. 

The integrand function X = X{Eo, Eq, s, E+, E_,^, rj, 0+) is expressible in the following 
(possibly the only one) viewable and compact form that facilitates the programming 
and the calculations to be performed (see also [|^, [1T6|): 

X = XiJ + Xy + X\Y^ 

where 

Xu = U + SiU + S2U + S3U + S2S1U + S3S1U + S3S2U + S3S2S1U, 

Xy = V + SiV + S2V + S3V + S2S1V + S3S1V + S3S2V + S3S2SiV, 

Xw = W + SiW + S2W + S3W + S2S1W + S^SiW + S3S2W + SSS2S1W; 

Si, S2, and S3 denote the substitution: 

Si = k^-k, p^^p;, pi ^ -p^, 

S2=Po'^ -P+, 
S3=p;^pl, 

where (with respect to the laboratory frame) k = {k,eo} is the four-vector of the 
incoming photon; p^ = {p^, Eq} is the four-vector of the incoming primary electron, and 
Pr = {Pr,Er} is the four-vcctor of the recoiling primary electron; p^ = {p^,E^} and 
jF = {pi, E_} are the four- vectors of the produced positron and electron, respectively; 
ko, Po, pt, p+ and pi are the corresponding three-component momentum vectors, and 
Eq, Eq, Er, and E_ are the corresponding energy values. The module of each 
three-component vector if is denoted by v. 
The expressions of U, V, and W are: 

U = i(l/(l + ny){il/kl)[-k3ikiT2 + koa^) + nk3 - r2(Ti - ra^a 
+ kiT2 + koas - k2k3 + k2 + n + 2k3 - (T2 + 2] 

+ {l/{k2k3))[a2{ki{T2 + ra) - 0-1X2 - 0-3X3) + k2{criT2 + a^T^ - 2x30-1) 
+ cr2{Ti - 0-2 + 2k2) - koki - Tik2 + 2o2 - k-^. 
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V = |(1/((1 + Ti)(l - a2))){{l/{koh))[2{ko -ks- 2ts) + ko{k, + n + a,-a2 + a^) 

+ k'i{—k2 - Ti + T2 + (72 - 0-3) + T3(-/Ci - A;2 + 2(72 " 2ti) 

+ A;o((Ji(-(72 - T2) - 2az{k^ + T3)) + k^{TiT2 + (7ir2 + 2(73X3) 

+ T3(2(t2(7i + (73T3) - kiT2 + /C20-l)]} 

- |(1/((1 + Ti)(l - (72))){(l/(A;o^2))[2(^o - ^2 - 2T2) 

+ A;o(^l + Ti + (73 - (72 + (7i) + k2{-k3 - Ti + T3 + (T2 - (Ti) 

+ T2{-ki -k3 + 2(72 - 2ti) + ko{a3{-a2 - ts) - 2(Ti(A;2 + T-2)) 
+ A;2(rir3 + (T3T3 + 2aiT2) + r2(2(r3(T3 + 0-1X2) - fciT;^ + k^as)]}, 

W = |(1/((1 + ri)((7i - l))){(2/A:2)[2A;iA;2r3 + k2{-ko -k, + ks-n+rs + a,) 

+ 2t3((T3 - ki) + ko + ki- 2k2 - /C3 - ti - T2 + T3 + (7i + (T2 - (73 - 2] 
+ (1/ (A;2A;o))[2((73(/c2T3 + A;3T2 - ko{a2 + T3) + 2x2X3) + T3(A;2Ti - /C1T2)) 
+ 2ko{ki - k3 + Ti + T2 - T3 - \ai - (72 + (73) + A;i(2A;3 - + (72) 

+ A;2(2t3 + (7i - 2ti - 2T2) + /C3(ti - (73) - 2t2(ti + T2 - T3 - (7i - (72 + (73) 

+ 2{ki - + ko-k2- 4t2] 

+ (1/(A;2A;3))[2((73(/C3T2 - A;3T3 - A;2T3 - koa2 + 2(72X3) + T3((7i/C2 - kia2)) 

+ ko{2ki + 2ks - (7i - (73) + A;i(2A;3 + T2 + r3) + A;2(ti - 2r3 - 2(7i - 2(72) 
+ A;3(ti + 2(t2 - r3 - (7i - (72 + 0-3)) 

+ 2(72(-Ti - T2 + T3 + (7i + (72 - (73) + (-2^0 - 2ki + k2 + k^ - 4(72)] 

+ (1/ {koks))[4:a3{k3 + r3)(A;o - rg) + 2r3(ri + ra - rg - o"! - C72 + (73) 
+ koi-2T2 + 2r3 - 2/C2 - 3(73) + A';i(2A';2 - + 2r3 + aa) 
+ A;2(2A;3 + n + 2t3 - ai) + A:3(-2r3 - 2c72 + 8(73) 
+ 2ki - 2k2 + Sks - 3ko + Ar^]}. 
They are functions of the invariant products: 

ko = p^.k = -s, Ti = p^.p; , (7i = p^.p; , 

ki^p;. .k, T2 = p^.pT, (72 = pi^T.pT, 

A;2 = P^-k, Ts = p^.p:^, (73 = p;.pl, 

ks = pqr.'fc, 

which are connected, because of the energy-momentum conservation laws, by the 
relations: 



(73 = fco + 1 - (7i - (72, 

a2 = ko - ki- Ti, 

(^1 = ko - k2- T2, 

Ti = /co - 1 - T-3 - r2, 
ki = ko - ks - k2] 
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where 

ko = —Eoeo{l — j3 cos 9) = —s (section 2.1), 
^2 = p-k[y cos 6k + (1 — y'^y^'^ cos(/)_ sin 6k] — E^eo, 
= p^k[x cos 6k + (1 — x^)^/^ cos 0+ sin 6k] — Ej^Eq, 
T2 = P-Po[ycos6o - (1 - y^)^/^ cos 0_ sin 610] - E^Eq, 
Ts = p^po[x cos6q — (1 — x^)"*^^^ COS 0+ sin 6*0] — E+Eq, 
(3 = po/Eo, 

cos 60 = COS 6 cos 6k + sin 6 sin 6k, 
sin ^0 = sin 6 cos 6k — cos 6 sin 6k, 
sm6k = {l-cos^6kY^^, 
cos6k = [k + po cos 6)/ Ptot, 

0- = 0+ - 

a=[{l + s- Etot{E+ + E_) + E+E^ + Ptot{p+x + p^y) 
-p+p^xy)/{p+p^{l - x^f'\l - y^Y/% 

Ptot = ipl + k^ + 2pok cos 6)'/^, 

Pk = So, Po = {El - l)^/^ 

P- = {El - if'\ p^ = {El - ly/'. 

With 6*0 and 6k we have, respectively, denoted the polar angles of the incoming 
electron and photon momenta with respect to the axis along Ptot = Po + k . In this 
case 6o + 6k = 6, where 6 is the angle between p^ and k , i.e. the angle of collision. The 
azimuthal angles of the produced electron and positron are 0_ and 0+, respectively . 
They are accounted for in such a way that 0^ = and (pQ = it, where (pk and 0o are the 
azimuthal angles of the incoming photon and electron, respectively 0. 

The variables ^ and 77 are introduced by the substitutions (section 2.1, equations 



(0) and (0)): 



X = X{^) = ^/l + Xmin (^ > 1), 

y = y{^, v) = bmax(a;) + y^U^W]/ (1 + v'^)- 
The remaining designations are given in section 2.1. 
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